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TITLE: Multiplicative signal correction method and apparatus 
Abstract Text (1) : 

An improved method and apparatus are disclosed for processing spectral data to 
remove undesired variations in such data and to removS interf ering information 
present in the data. The method land apparatus corrects multiplicative effects 
present m the spectral data. Additive and interferent contributions can be 
corrected as well. In one aspect of the method, coefficients for a selected 
appropriate model are applied to the input spectral data based on first and second 
reference spectra. The spectral data are then corrected based on the estimated 
coefficients at least as to multiplicative errors for producing a linear additive 
structure for use in calibration, validation and determination by linear 
multivariate analysis . The method and apparatus will iitprove the accuracy of 
spectral data structures derived from measurements Using spectroscopy 
chromatography, thermal analysis, mechanical vibration and acoustic' analvsi" 
rheology, electrophoresis, image analysis and other analytical technologies — '~ 
producing data of similar multivariate nature. 

Brief Summary Text (2) : 

The present invention relates generally to processing of spectral data to reduce 
undesired variations and to remove interfering information present in the data 
specifically, the present invention relates to an improved instrument or method for 
processing of spectral data to reduce undesired variations and to remove in^erferina 
information present in the data. Most specifically, the present invention re lates to ' " 
an improved instrument, method or process to provide improved measurements of 
analytes based on spectral data by reducing Undesired variations and removing 
interfering information present in the data. 

Brief Summary Text (4) : 

Spectral data consists of multiple interrelated data points, such as an optical 
spectrum or a chromatogram, which carries information related to the components and 
characteristics of the specimen from which the data was derived, as well as to 
performance of the analytical instrument and to the general experimental conditions. 
In spectroscopy, for example, this specimen is a material and the spectral data 
comprises the results of related measurements made on the specimen as a function of 
a variable, such as the frequency or wavelength of the energy used for measurement 
In chromatography, the spectral variable may be time or distance. In thermal 
analysis, the variable is usually temperature or time. In mechanical 
vibration /acoustics analysis the variable is usually frequency. In rheology the 
variable can be position, shear rate or time. In electrophoresis and thin layer 
chromatography the variable is relative distance in one or two dimensions. In many 

.^-g- kinetic measurements, time is either the primary variable 
or an additional variable that adds to the dimensionality of the data. 

Brief Summary Text (5) : 

In image analysi^^^ the fundamental variable is usually distance in one or two 
dimensions although the two-dimensional Fourier transform, also known as the Weiner 
transform, and the Weiner spectrum which express the information in a 
two-dimensional spatial frequency domain are also prevalent. Multivariate images 
such as three color video signals and many satellite images where each picture ' 
element is characterized by a multichannel "spectrum" and also images constituting a 
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AitLSeS^ln'iSr'^j^^^^ dimensionality in the data, 

showine dist?lbutioro% 5iH^ic\teTSen?^ 

^lliliiMilF"- 

.^ses ,1,0 apply to th,ee-way to.o3„p.ic i^sea.l'a' Jr^S'Tna °x' ;aT"t'S;^raphy. 

Brief Summary Text (7) - 

iTrllltll t'o%:crotJe'?%%a°ti:?Sni to I: ^"^^^-^^ -^eral specimens .ay 

derive specific ^e^lT/ .11^^^^^^^^^^ coSSi??on'or'^^ '° 

constituent, the magnitude of a physical- property, or the JdSntmSio^^^^ 
specimen or its components, one form or another of aSdiMve ™ufia?ia?e 

Brief Summary Text (8) • 

spectrum, or estimated "loadinq" vectors (e a i i r^c.^Z t ^Zt 

spectral effects, thereby providing selectiJitJenSanceLn? ' ^"terferants ■ 
Brief Summary Text (9) • 

rby'L\rng'tJ:^'Sri?K'^Se°Sat; ^T°^"f ^ ^"-1^ ^^^^itive structures, 

»lupucauve .,o„„, pai Jcu^la'ir^^H^^ISS Sa^-^^naKTlf 

Brief Summary Text (10) • 

Near-infrared spectra ^enS to dec^eL?Tn^h. k ^""^^ examination of spectral data. 
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order of the harmonic involved increases i p i-v,« e 

multiplicative effect on such a til^^H oAo^; frequency increases. A 

tilted baseline There?oreSorrisJS o?hf ^^^^^^^ ^""'^^^ ^° addition of a 

Brief Summary Text (11) : 

Removal of a multiplicative factor imnlies diviq-inn of t-y.^ w 

value. Norris (K. H. Norris and P r wtn^f^o n ^ ^ "^^^^ appropriate 

H. Marten, and H. Ru,=„™, Ed. AppUed SleS^ki^rf L?d lies S=L ^°^>'"^ ^ 

"."°L5''.L'L'i„r„??:aL'?s'iJL^s.°?nr,rnS'^rs^™iSr^^^ 

parameter of Interest . A dJZi°ItS\fve!e„"th'ls th» Seot^d'Sv'sJJflr °' 

dii~~^ 

difference computation of the derivative is then performed to oD^imiL thP 
approximation. Additional terms may then be added ^o hh= ° ^ ® 
procedure. This method has been ^Jeful SowevJr Jt is Umiteflr.", ^ "'."^"'"^ 
calibration using data at a few selected wlve^;ng^hs ° ^ '^'"''^ 

Brief Summary Text (15) • 

aSoJelor^rof speSrlJ ZtT'^^r.r.ft'''^' interferences also characterize the 

specimen or instrumentation, and stray liqht sdscuI^t- r-»fio^^?^I^ ^ 

SS;^eed'l^^'rL°"',JL^af;L..e. 

Brief Summary Text (16) • 

speJtral'^f^ec^s'of'sSci^aria^bJL^^iitL^f'list^^ "^^^ 

to use of the spectra for calibration v^H . ^ analyte information prior 

the utility and^perfo^°n^crif^L1^ie;rnie'S:"' SlSlTech^er^' '^""'^^ 

Brief Summary Text (17) • 
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Brief Sunmnarv Text (20) ■ 

spectroscopy, "^hSil^rLhv SeL^f ^^^^i^ed from measurements using 

anal^ '^.eolo^r e?ec?ro^horesT^ ' i:.m^.:n''^lTo.M^^ -^"^ 
technologies producing data of simil'arStiviHali^ture ' ^"^^yt^^^^l 

Brief Summary Text (21) - 

rl<ltle'or°lir.^^^^^^^^ correct spectral data to 

subsequent additive modeling ^"^^ improving and simplifying 

Brief Summary Text (22) * 

multiplicative signal correction process . clesired information m the 

Brief Summary Text (24) • 

interactive „o„-U„ea. aodelin. i'L'^:!".%l7c.f.'J^AlZr^^^^^^^ " 

Brief Summary Text (25) • 

Brief Summary Text (27) • 

iorrLSn^S'^SS^'^yeiJirStel^LtS^r ^'^^^^^^ ^"'^ multiplicative 
interactive as well as fuTlv lula^tt^ nn . Provide graphically based 

their judgement and ll.lr^iLl^'lr^^^^^^^ - ti^y^sf desire" ^° ^ 

Brief Summary Text (29) • 

correction method in t-hf* f-ive^ ' ^ - prior multiplicative scatter 

Brief Summary Text (31) ♦ 

steepest aescent »ethcl ^™art?.; n^.^ ° ''"I""' ^.tpansion, the 

coefficients ^ni^i^i^'tre'^^ S^r^r^J the^Zl'^e!?""'^"" ^° "^"^ 

Brief Summary Text (32) • 

'^^lnTl^'e%7ss'::r^^^^^^ --^^^-^ the multivariate 

BFthrlnterf exing ^' "^^"3 the coefficients 

the spectra • of Thls^ c^oneSs to St sub^r.ot"^^'' 7 non-linear rnodeUng to scale 
data can substantially ^S^vrt^SL'^oStSL^^S^fJof tL'^da't:'^' ^^^"^^ '^^^ 

Brief Summary Text (33) • 

ip^?^nrtirss^s;°cSiK^^ 
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oni^i™''^^®T''^ ^^^'^^''f ^^^"^ ^'^^ orthogonal to, and therefore uncorrelated with 
one or more reference analyte spectra. The coefficients generated for ^heSe 
?S™Mi """"l"" are not influenced by the presence or^magni^ude of aSlyte 
information contained in the raw data even if the analyte spectrum is not included 
in the modelinc) or the coefficient estimator does not inhe^entlTorthc^onaUze^^^^^^ 
JpSr^J'^ff^cts ' t representation of the LgnU^de of the'"' 

hh!n ,^=L ; , u ® ^"=-^^gerinq components. These more accurate coefficients are 
then used to scale the original reference spectra prior to subtraction frnm^h» 
input data and the correction proceeds as in option A TMs op^Jon reSuces or 

^"^'^'^ spectral contribution that woSld otSerSJse be caused 

iii::^.r.Tio\:^^^^^^^^ °^ ^ --^ -tains s^Tniit^tiir^' 

Brief Summary Text (34) : 

I^'I P^^^^^^ i^^^rition also comprises, as option C, the further scaling of the 
spectra of the intertering components and, if desired, the spectra of the analvt^fO 
data' Soec^^a^' °' modification aAd correction Applied to tJe^'^ 

br?onSrS\rSe°SEt\^^^^^^ down-weighted, left as is. S e^hasized 

Brief Summary Text (35) : 

The present invention also comprises, as option D, updating of the analvte and 
incerference spectra based on the results of later s'^ages of data So'essTng and 
^^^-y^^s for example based on principal components analysis (PCA) (S^'nartTal 

*h\"- particularly useful in conjunctlS^T^ ..S^^nal "^^o^L.T 

Method and Apparatus" U.S. application Ser. No. 07/319,450 filed mI^ 3 I 'isThl 
Edward Stark, one of the co-inventors of the present invention 

Brief Summary Text (36) : 

oI!^r,f^^^"'' invention also comprises, as option E, interactive displaying graphical 
output concerning which analyte and interference spectra, if any aL cauJtnq 
difficulties with respect to estimation of the multiplicative corJ;ctron and ^ 
interactive control over which reference spectra are utilized, the sp^c^ral ranae 
^™ons ItlTylr weighting of' the ^^"^^ 

Drawing Description Text (5) : 

FIG. 4 is a block diagram of a non- linear model coefficient estimator used in 
conjunction with the present invention. 

Detailed Description Text (2) : 

The present invention has general applicability in the field of signal and data 
processinq, wherever -spectra" or data structures consisting of multiple 
d^^r^'h f^'' data points are obtained and the variability in the data can be 
Trl rnllt ""'^^"f additive and multiplicative effects. These types o? effects 

Detailed Description Text (4) : 

Like Murray and Jessiman's "mathematical ball milling" it seeks to normalize everv 
input spectrum to some reference or "average" state by addi?Jve and m^T^ionc^^r^e 
normalization and it allows an explicit correction for wavelengJSef^ecJs (mosJ 

efS^iveli"ti™S ^Sr''"^'J " '"'^^^ "^'^^^^^ "constituen?" vCctor but Sore 
effectively through the use of non-linear modelling ) . m addition, it employs a 
different regressor (an actual reference spectrum, e.g. an average ^pectS for 
determining the multiplicative correction and it allows nrodeUno and ?f dSr-^H 
subtraction of several additional phenomena at the same ti, !!^. ^ ' 

Detailed Description Text (5) : 

\i!S:)^^tt''tZtTtT^^ Martens, Jensen, and Geladi multiplicative scatter correction 
"av^J-^nt" It^ I normalize every input spectra to some reference "ideal" or 

T ^^<3itive and multiplicative normalization, after having estimated 
=nL^ ^""^ parameters by some type of regression against some SfeJeS^ 

Ts lZ"^ w^^H ^elected wavelength range, andlhis reference spect^am S^Se of 

the same kinds as those employed in MSG. However, it extends conventio^ST S^C by 
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. ■ b«P*»al..80O2W8a..«.?im^ .,_2^ ,t«j^ to 4^ ik» fat, fn, 

explicitly «0d«llnq the effects of anticipated additive in, „f and hv 

It allows removal of nrilTi^lTl^^Z^ SelL ^t tZ St'^"'' 
?oSect-raL1-Jli?a^ri?^"TSS4^^ 

contribution to ThS^fSStSl data b"r^5"„:Ji;i°,2 "^''^'^ -"'"ction of their 

Detailed Description Text (6) * 

Detailed Description Text (9) : 

This physical situation described above is linear in the nar;,m«^ar•=, • 

a^Klu" tfoS: Ttifva^LS: r°X-l "^ l^r.^^L that 

that is a function of k can be consideJed an fdSt^fon i ^^'^ f"^ additive baseline 
Variations in the additive of f set aJd tS ^ /^f'^''^^^ component Q. sub. km. 

from different specimens contributes the adSt i^anJ mult^oirc^r"'"'"' """^^ 
which this invention provides improved daSa JSectJons '^''"^ 

Detailed Description Text (10) ♦ 

rsliTF"-'.'?--™^ 

the input data that is not represented by any of the ?eSren=e sSSa Jhe''"" 

S'-e"tS:tTr.r.^^-;b-^^ 

Detailed Description Text (13) • 

th^analyt:: "rinteSt on^'tJe'^L^u^rSata "^^^-^^^^"9 --P-'^ed influence of 

to indicate the .u.'llly s^^s^.Vnl^l su^sec^en '^^naJvS TSe' 

?SCien^tf ba%r:^^-b^- — ^ on^ ?r' mtVoTTt^s^ ^ 

by the instrumentatiSir'^rihe'il^E^^sib.ki' °" artifacts introduced 

Detailed Description Text (15) • 
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characteristics of R sub ks be stablP if I k ^^^^^^^^^^ that the spectral 

solvent spectrum are iSluded S r sub kj ^ devotions or variations In the 

Detailed Description Text (16) • 

orthogonal sets of spectra re preslntlt ive of t-hl IX . ^ squares (PLS) provide 

span the variation of X.sub.ki so as to stabilize thr Li!i reasonably 

and specificity are not essential For ^h^^o f "'"^^^^'"'^ ^""^ spectral accuracy 

Detailed Description Text (20) • 

iiiSiiaii* 

in transmission spectroscopy, where an increase in rpHnr^=;h convergence error 

Detailed Description Text (21) • 

rS-Se\j^nn-Tl£!^-^eL"^S"rt-L^L^^^^^^ 

Detailed Description Text (25) ♦ 

Detailed Description Text (26) • 

The corrected spectrum Y.sub.ki comprises the standard spectrum R.sub.ks and linear 
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additive deviations from R.sub.ks caused by analytes, interfprants and errors »nH 
It IS therefore suitable for further linear data analy sis . ^"'^ errors, and 

Detailed Description Text (27) : 

Preferred embodiments of the above are illustrated in the figures and further 
described below. In FIG. 1 for instance a photospectrophotoml^rt; sensor system 100 
as used m the present invention is described. This system can be used forSmoJp 

fiber optic bundle 130, which transmits it to the spectrophotometer 140 ?hir 

fn^nf'' ^""^^"^^ '^'^ 1^^' a concave holograpMc ' gracing 142 

and one or more diode array detectors 143 and their associated o?de? sor?S JflJer; 
144, arranged to measure energy at different points in the spectral iS formed b^ 
ana ^ea^^inS.'" f ^"J '^Jj-^f^^^ ^t different wavelengths with^the v^sSle"^ 
and near-mfrared regions of the electromagnetic spectrum. Each detector channel ha^ 
an associated preamplifier 145, the output of which is multiplexefby ^ItSexer 
furth2'°H^ programmable gain and offset amplifier 147. The spectrophotome^e? 1^0 is 
further described m "Improved Grating Spectrometer", a U.S. patent apDlicltion 
filed Aug. 24, 1989 by Edward Stark, one of the coin;entors of the present 
invention. The application has issued on Mar. 5, 1991 U.S. Pat No Tgg? 281 tHp 
contents of said patent are incorporated herein by reference 4,997,281. The 

Detailed Description Text (28) : 

nLtl fT^ ""^^ ''^"'^ sequential multiplexed analog signal is then 

sSst^m T ^" to digital converter loi in dati ac^isition 

system 200. In preprocessor 202, the energy data is proc^-oed to elimin^Tt 
instrumental offsets and to reduce both systematic and^lsS^oise and the^ 
to obtain data relating to transmission of the specimen. This data is then 
linearized with respect to the analyte information of interest eg Se logarithm 
of transmission IS more or less linear with concentration of chemiLi consti^ue^^s 
th?i^^^ H ^Pf^;"^"-.'^*^^ data is then in form to be normalized in accordance 
the methods of this invention. Although the details may differ, similar functions 
are utilized in obtaining spectral data of the other forms discussed above 

Detailed Description Text (29) : 

The additive and multiplicative corrections of this invention are performed bv data 
n™ "^t-'^^'.f comprises special purpose digital computation Joglc After 
normalization, the data may be further processed in processor 400 prior to use for 
Zl^Z '^''''^^'T' ^^li'^^ti^" and d etermination^f unknow n values For 

example, processor 400 may comprise the invention of "Signal ProrPssira Method ;,nH 
TontT.T\''-^-/^^T ^PPli-^tion Ser. No. 07/319,450 filed Mar. 3', 1989" Thl ^ 
contents of said application are incoporated herein by reference, plnally the data 

ca ugSfon v*"; S^""^ "^^^'^ functions of ailtivariate 

calibration, validation, and determination required to generate the analytical 
values then presented on display 600. analytical 

Detailed Description Text (30) : 

The data normalizer 300 of the present system illustrated in PIG 3 provides a 
number of options for processing the input spectral data X.sub ki The basic 
orre?erencf°'' °l ^°^^'--^-on provided by this invention is based on the use 

sJecS s?oL5'? ^t^'^d.^'^th'^'' ""'^ ^■^"^•'^^ ^"^^-^d - the reference 

^^I^^ • ? ^ . ° ^"P"*^ spectral data X.sub.ki by means of the 

coefficient estimator 320, and to determine corrected spectral dat^ J sub ki bv 
wht:!h ^^l<^"l^t°^330. Functions 340, 350, and 360 provide adStional options 
which are bypassed for the basic corrections. The control and logic sequencer 370 

si?ecf^H f ^ ""^'^ selection, and control signals required to pS?om the 

selected functions in proper sequence. = m i-u pertorm cne 

Detailed Description Text (31) : 

i^i^ preferred embodiment shown in PIG. 3, this correction is implemented bv the 
calculator function 330 comprising subtracter 331, divider 332,^Gbtractor 333 and 
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divider 334 that perform successive operations on X.sub.ki involvinq the 

Si^t^r^lo'^lhe,"'' aeneS Te coefficient 

estimator 320 These operations are performed sequentially element by element bv 
indexing k with a first counter and performing the required se^ence of drq^taJ 

?hLf 5igitararj?SLSf f 'T"" ^^^'^^^^ '^^^"^ °" ^^^^^ ""^'second c^^^t^r. 
inese digital arithmetic functions are available digital logic functions utilized in 

an^raTi>S^:^°"^'"'^""^ ^'^^^"^'^ ^" n^h^coprociss^diiice^'or" 
Detailed Description Text (32) • 

In ?ir1^''^HHif'''^''?°" '° """^ ^^"'^ preferred embodiment described above and shown 
in FIG 3 additional correction spectra [R.sub.ka *b.sub.ai ] and [R sub ki 
*b.suboi ] are formed by matrix multiplier 340 from the reference spectra and th.ir- 
associated coefficients that were also generated in coefficient estiS^tor32f T 
matrix multiplier to form the spectrum [R.sub.kn *b.sub.ni ? Sr f sS?e inout 
spectrum i, consists of short term storage for the both innuts ^ n,nt^?r;?- T 
sun^ation Circuit, an address sequencer whicrac^e^ses^i^'coJ^e^pSiif eSentsT 
of R^sub.kn and B.sub.nk and a second address sequencer which accesses ?he rowr^ of 
R sub.kn and addresses the short term storage wh?ch keeps tSe resulting J t^mes 1 
matrix. Matrix multiplication is also a standard function of aJaUaMe a^rat 

^^^2l^'LT\T''"'- correction spectrum may be used direcUy by 

subtractor 333 to further correct Y.sub.ki, which becomes 

Detailed Description Text (34) ■ 

artCis°siLe °^ """^ possible [R.sub.kj *b.sub.ji ] corrections are applied 

eff^Mvf th^n ^hf Processing, because later additive linear modeling ^y be more 
effective than the spectral subtraction based on previously known referenc4 soectra 
performed here The [R.sub.kj *b.sub.ji ] should include, however, those 
^":-^^^^""ts that are difficult or impossible to adequately represent Tn the 
calibration data, e.g. moisture and temperature as previously discussed 

Detailed Description Text (35) : 

Greater control of the situation is provided in a preferred embodiment that also 
incorporates multiplier 351 and component weight storage 352. In tMs eSodimen? 
foLfno^'f °K r'^*'" °" °' "^"^ correction spectrum that is s^btrStSd in 
obSJned ir^ater'data' ^^^^^^^^^^ ^P-^*^- - accordance with i^ioJILtion 

°b sirfi ? in^ h-i^p Processing steps. The corrections then become ha[R.sub.ka 
MElmM h3 [R.sub.kj •b.sub.ji ]. The corrected spectrum then becomes 

Detailed Description Text (36) : 

J/n^^?^ °5 ""^^^ "^^^ complete cancellation is desired, while a weight of 

h:=0 provides no correction for interferant j and ha=0 prese^es the analJte Jianal 
unchanged Values 0<hj<l may be used to downweight infornStloHhlt hL ^ncfrLiX 
or potentially harmful effects on later data analysis without toLfJeiicSn ha <0 
increases the weight of the analyte informatio n, thereby reducing the SSi? 
importance of other information in the corrected signal relative 

Detailed Description Text (37) • 

^"^^y^^^' they are orthogonal by definition, if they are measured^ectra of 

components or otherwise separately derived, they will generally be fnter^^^^^^ 

which If severe enough may cause errors in the Coefficient valLs or failure of th^ 

Detailed Description Text (45) • 

Is\orJ?fiTH°K ^•^"b '^\^r« preserved in P.sub.ks, therefore the coefficient bs 

^"ected by the orthogonalization. Each succeeding R.sub.kn is theJ 
orthogonalized against the matrix formed by the precedini orthogoLi P suS kn 



9 of 20 



10/23/2003 1:39 PM 



Record Display Form 



importance, or to accept it without change OrthLonSfrir.- on '^""'"^^^S^t its 
performed solely for the puroose of aen^r^^ino processing may be 

orthogonalization is chosen^the ^eferenJe spJct^a innJ^S ^f^^^^?"- '^^^ 
are the P.sub.kj and P.sub.ka. ^'^^^"^e spectra input to matrix multiplier 340 

Detailed Description Text (46) : 

The orthogonal component generator and storage 360 comprises stor;,a^ fo^ d k 
the portion that is filled a<3 tha -^v-,-^«„„ -^o" >-umprases storage for P.sub.kn, 
rz-71 1 c^-J^^Z I ^^^^e<^ as the process proceeds comprising Z, storaae for 
IZ ZJ-1, storage for the intermediate Drodur^ 7i7<7\ Z,^ Z: '*"^°^^9e ror 

for the intermediate product Z • Z subi no^n^ ^'^^ storage for Z.sub.i, storage 
scalar inversion (1/a) Jo^ic a suStra^tor anH''^^"'' "^^^^Pl^ sum logic, 
storage for process' ng to Control ^hf " ^ sequencer to select data from 

devices (PAL'S) coirtoined with a M^iu ^' standard programmable array logic 

sequential fashion by the se^enJer. ^^^^^ operated in 

Detailed Description Text (52) • 

5. Move spectrum R.sub.kn to Z.sub.i storage (R.sub.ks for n=2) (K.times.l) 

Detailed Description Text (64) • 

d^fficiJJ^^r^ejSnj:: ^liTtiZiTlnTtl'''' ^° drastically that it becomes 
aliased compared ?o the'S^Jginl '^uaSt r^p^^s^nted'^bf the^'^"?^ thoroughly 
This is particularly troublesome 2hen inte^^^nce 's^^^^ 

downweighting, or analyte enhancement is-^^I^T^esf Stnr^Af ^''''^"^^"^"^^ ^ 
desireable to perform less drastic processing factors often make it 

Detailed Description Text (65) : 

An alternative preferred embodiment orthogonal i 2 es earh -r-ov-f^ 
against the analyte spectrum by the simplfuj^^rrgressioir^l^ 

Detailed Description Text (66) • 

described previlusjy ^ requires only a subset of the functions 

Detailed Description Text (68) ♦ 

t2\£r"Jjo:fc::pj:^;S used":o'Todtr\t'^'T'H"^^ -efficients applicable 
estimation orocasa involves creatin'^^^^^^^^^^^^ '^^^^^ the coefficient 

data that is^rfunction of theTeferen^^ iS'traf da^fanS'"' °' '"^"'^ ^P^'^'^^^l 
other variables such as the input spectra^'^datritsel? Sd'k!" "^^^-^^^ lOe^^ls,. of 

Detailed Description Text (69) • 

L"cS1o"r:pre"S o'si^^^i^ :'lkluT^ lVlr'' '^-''r''-''i * 

coiu™ is , reference spectruH . ^^^^J^L^pSL^e^d T'^^ !Sr'ci"lf i^itt 
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estimator 320 fits X.sub.ki to Rnk by some method, minimizing the residuals e.sub.ki 
Detailed Description Text (70) : 

Methods for achieving this linear cnodeling include generalized least scniares 
maximum Ukelihood regression, robust regression, estimated best llniar Tedictor 
partial least squares, principal component regression. Fourier regressLn 
covarxance adjustment, and others. For example, gener;iized least Safes' with 
generalized inverse models X.sub.ki by squares with 

Detailed Description Text (72) : 

However, a preferred embodiment uses the more usual linear modeling performed bv 
multiple linear regression where — periormea py 

Detailed Description Text (74) : 

can be precomputed externally and stored with the reference spectra, thereby 
minimizing the requirements on the data normalizer 300. If full rank Gram-Schmidt 
orthogonalization is used, U.sub.nk is available from that In either case 

mulMnV f the coefficients of the linear model invSwiTT simple matrix 

multiplication. A matrix multiplier for U.sub.nk and^TiUb.ki consists of sho^t term 
HZT. °^ '^^'^ '"P^'^"' ^ "^Itiplication and summation circuit, an aSes! 

sequencer which accesses the corresponding elements k of U.sub.nk and X sub ki and a 
second address sequencer which accesses the rows n of U.sub.nk and addresses the 
short term storage which keeps the resulting b.sub.ni values In the moSqene^ll 
obt'Ln mVif^P'^ ^J"^^^ regression, matrix [R-R] must be formed aS in^^rted to 
obtain [R'R] .sup.-l prior to matrix multiplication by R- to obtain u sub nk This 
iSJc'sequ'e^ceT'''" accomplished with an availabL array pro "essor 1nd'su?tah1 . 

Detailed Description Text (76) : 

oL^?! "^^^ of nonlinear fmodeling, the coefficient estimator 320 becomes more 
complex as each nonlinear coefficient becomes a vector of length k C sub ki and 
D.sub ki are therefore matrices containing a number of coeflJ^ient 'vectors Jha? 
depends on the form of the nonlinear model . vectors that 

Detailed Description Text (77) .- 

These coefficients can not be determined by multiple linear regression or other 
bilinear methods so an interactive procedure must be used. Methods in the literature 
include linearization by Taylor series, steepest descent Marguardt • s compromise and 
simplex optimization: (N. Draper and H. Smith, Applied Regression Ai^aTsTs Second 
Edition, John Wiley & Sons, New York 1981 pp. 458-465) . ^ ' i-econd 

Detailed Description Text (78) : 

A preferred embodiment uses the coefficient estimator 320 illustrated in FIG 4 
which employs Taylor series linearization. The model response generator 321 
calculates the vector F from the reference spectrin, the presfnt Sue Ari of the 

spLt^^rSata'^xj;^ 'ZTolTr t^ 'T^'"''' P^o-^^^^- the'variable k and th°e inp^t 
«™^=r;^ operation involves matrix multiplication and summation in 

accordance with the appropriate form of model as discussed above. The^t of 
coefficients Ari, comprising cni of Cki, dKrEf Dki, and bni. aSe initiaUy stored 
JSrSh aSd!?ioJ'';f ^^°^^9%^22a. They may be modified by me^ns of add" 322b 
through addition of a weighted correction wGq or of increment dAr to one of the 

^.^^^''h^"^'/^ ^ "^^"^^ P^^^^"*^ -^1"^- stored in coefficient Aq storaoe 

obvLus f^om tS^ ^.response generator 321. The remaining functions wi?? become 
obvious from the following description of the operation. 

Detailed Description Text (81) : 

3. Regress Xki on Rkn to determine the linear model bni; 
Detailed Description Text (93) : 

iLi^^'i'* obvious to those of ordinary skill in the art that some of the above 
operations may be performed in different order without significantly affectina ^h» 
results Obtained. It should also be obvious that the func?Jons sS cfn hS) f be 
implemented with common digital logic circuits well understSd by S^ose of ordlLr^ 
Skill or by available microcode controlled array processors; such as the Data 
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Translation Model DT7020 with the MACH d^p c;„hy.^„*-^ „^ t k ^ . 

of the appli^iblTData Translation ^Znss (S'u ss ^at^^c^^^^ ""T^^ 
pages have been annexed hereto and inco^orated " li '^^^^^^^^^^"^e^^^ScT 

Detailed Description Text (94) • 

irtrafaradlJi^''^ method used to estimate and correct for multiplicative effects 
fasMon ^"'^ ^ multiplicative in an interactive sequential 

Detailed Description Text (95) : 

1 Let X=(x.sub.ik) be the matrix of spectral ordinates for i=l,2, . . n obiects 
Sa1;^;sin, a st\S^S JplL^rj^^^SS^™^^^ from tht^^p^e^c^^.al 

components for practicing tSe present ^^^^^^t^^nTSJd^n^ ^he ' c^^rec^^ s^^J^ll^^ 

Detailed Description Text (96) : 
2. Estimate an additive model : 

Detailed Description Text (97) • 

P^!p%nh'^f f °^ ^^^^ ^' ^-s^b.mean is the mean vector of Z and 

P- (P. sub. kl) spans the spectral variations of analytes and interferences a^ w^T i 
possible and/or to the extent the user wishes. P may include a^T^^^^^i^- 
input component spectra, estimated component spLS/JoaSJngrLom f PcJ°or 
fS^LS v"" ''"'"3 estimated component spectra or loadings f rom 

l?..ZefV-l^,^l i°r%Lm?ir"^^"^^ ^^^^ °^ -^^'^ 

Detailed Description Text (102) : 

In this method, the multiplicative effects, say from a physical model and the 

W^^^^^^^^B^BB ^^^^ 

Detailed Description Text (103) • 

b«J"of "»P"=''=^."='=>'Ol of data normalization provided by this invention 13 

based on the use of previously obtained analyte and reference soictra rj, 

^i''?""??"^ °" 'f'""' us, o? S iSSSn does 

<:r»^i— ?e^, 

=d ^tS™ ^t'?hTs"LS^?^~^ 

Detailed Description Text (105) • 

where X-U.sub^ik) is the matrix of spectral ordinates for 1-1,2, . . . n obiects 
Mnn;ar £act»s"!ST"°' '<» obJe«s H 

componn;t' °analvsl"- 'oa^;t.a, leistti^^, e't" r°^-r^"u"b"l, a^e 'r=d-^""f 

: rSei °T%^"rLid?„°7.''caT£'brd"'" 

SDectra r-fr «nh ^® decomposed into a function of a reference 

Detailed Description Text (106) • 

rengVK''c-(^sib'?m)'ri'a':i.'^ "^""'"'^ °' '^"^^^ ^' ^ - ^ °f ones of 

whV^^ L-ic.sub Im) is a matrix of regression coefficients of size L times M 
which quantifies the analyte and interference contributions and P-ff '^nh 
contains the residual loadings with WSiTt^^^^^^ 
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phenomena removed, d. h. and C can be estimated by regression of P- on r i ;,nH r 

Detailed Description Text (110) • 
Detailed Description Text (111) : 

C.sup^* and G.sup.* are user- chosen subsets of C and G which include tho^P 
xnterferents andanalytes of interest which it is desirTd to eliminate from ^h« 

Detailed Description Text (112) • 

directly to the multivariate scores rather than to the input spectral data 
Detailed Description Text (113) • 

Using^the prior example where the input spectral data X are modeled using a bilinear 
Detailed Description Text (114) • 

l%°lli%r:ntll^ ^nterferen^-corrected spectral data can be defined as 

Detailed Description Text (115) • 

analyte and_ interference information contained in C. In practice i? additive 
Detailed Description Text (118) • 

The product T*U produces offset and interferent corrected scores and V- i <, ^h« 
matrix of corresponding spectra loadi ngs associat ed with the corrected Lores 

Detailed Description Text (119) • 

Tc^n^lnT^llolTo^lll^'' iSEKlSSL corrected dat. z can be 

Detailed Description Text (122) • 

Detailed Description Text (125) • 
Detailed Description Text (126) : 

2. Using a statistical method (e.q. weiqhted Ipac,^ omi^.r-oo^ ^ ^ 

and appropriate analytes and in^^e^ G.'^c^cu^Tt^^^^^^^ ^' 
Detailed Description Text (127) : 

3. If the spectral scores are to be corrected calcnla^p n (i=r^^ ^a^-^- 
iLll-ferent. effects) and S.sup.-l (for nSitfpli?at?ve effects) ; 

Detailed Description Text (130) • 
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ilchSiS^s^^^' ^''^""P^^ component regression, Fourier regression, and other 
Detailed Description Text (131) : 

aDDrooria^^m^^hof ?-^f''r''K^'^ spectral scores W to fit a linear o-uodel using an 

appropriate method listed above. 

Detailed Description Text (136) : 

3. use the corrected data and the calibration model coefficients from the linear 
'"oae^ to predict the properties of interest. linear 

Detailed Description Text (137) : 

M aSf "^^f ^iP^i°" °f the alternative embodiments described immediately above (pp 
35-40), the apparatus described generally in PIG. 2 still is applicable as would be 
b?oc^f °S ''^ ^on-'^-^tion of some of the more detailed 

blocks, the coefficient estimator 320 described above is preferably the basic 
element. For example, the estimation of an additive model (p 35, step 2) is 

coefficient estimator 320. The reconstructi^JTof the spectral data (p. 
I' Preferably performed by calculator 330. The iteration required on 
«^ = ^o^ ^^^^ controlled by a logic sequencer 370 or equivalent. Modpling from 

statistical analyses (principal component analysis or partial least sq uares for 
example) may be accomplished by the structure shown in F^G. 5. DecomposS^Jn ^f 
loadings (see page 37) may be performed by the coefficient estimator 320 Other 
functions are readily performed by apparatus disclosed herein. 

Other Reference Publication (2) : 

Pood Research and Data Analysis edited by Harald Martens Hellmut Russwurn, Jr ■ 
Applied Science Publishers, Ltd 1983. uo»wuxii, ur., 

CLAIMS : 

Lr-t Zf^'^ ^°rr«'=^i'^9 Signals representing input spectral data (X.sub.ki) 

«fi^ .J''^'" ^ '^"f^"^ ^ measurement, at least as to multiplicative errors, 

said method cortprising the steps of: 

providing a first and primary reference spectrum signal (P.sub.ko) representing a 
predetermined standard for such data; representing a 

providing at least one second reference spectrum signal (P.sub.ka or P.sub.kj); 

f"*^^ f°r/ selected appropriate model to be applied to said input 
data based on said first and second reference spectra~iiiI7als; and 

correcting said signals representing said spectral data based on said estimated 
coefficients at least as to multiplicative errors for producing signals representing 
a linear additive structure for use in calibration, validation and detJn^iSI^ioSK 
linear multivariate analysis . ucucj-miudcion oy 

l\Jr!l^ "^^^^"^ ^ wherein said at least one second reference spectrum 

signal represents expected influence of undesired interferences (Pjk) ^ 

4 The method as in claim 1 wherein said at least one second reference spectrum 
ipTl IT^^T^ ■T^^'^l^t i"fl"^"<^e °f analytes of interest on the inpu? data 
(P.sub.ka) and said method also includes the step of providing at least one third 
In^lZll '•''f (P.sub.kj) representing the expectel influence of ^ar^ous 
~r'ie^&fm^ ^"'""''"^ '''' '^^'^^ "^^^'^ °" ^^'^ - °- 

tecLr^e"'* °^ ^ "^^^^^^ ^^^'^ includes a generalized least squares 

7. The method as in claim 1 wherein said model is a linear inode.! . 

tectargre^^°^ °^ ^ ^^^'^ ^"^1"^^^ ^ "^^i"""" likelihood regression 
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predJct^f SchnL'^^''" ' ^^'^ ^ -'^l"^- - --^i-ted best linear 

regresLor^SchnL'e^'"' ' ^^^'^ SE^el includes a principal component 

llcZT^T'''^ °' ' "^''^ ^^odel includes a covariance adjustment 

12. The method of claim 1 wherein said model is a non-linear model . 

13. The method of claim 12 wherein said model includes a Taylor expansion technique. 

14. The method of claim 12 wherein said model includes a steepest descent method. 
llchlT^T'''^ °' "^'^ ^"-^^"^^^ ^ Marquardfs compromise 



16. The method of claim 12 wherein said 
technique . 



includes a simplex optimization 



coefflciZl^ f \^ ' °^ ^^^^"^ ^ °^ ^^^i"" ^ including the steps of using the 

the sn^r^^r.? ^^'"^ ^"terfering component derived by the modelir.g ?o scale 

the spectra of the component and subtracting the scaled spectra fro m the data to 
remove their contribution from the data. f ^ <i trom tne aata to 

i^^J^l^ T^^^^i ^l^i-" 1^ including the steps of generating modified reference 

spectra of mteriering components that contain only those portions of orSinjr 
reference spectra of the interfering components that are orthogonal to and 
therefore uncorrelated with, at least one reference analyte spectrum, us'inq the 
coefficients generated from said orthogonal reference spectrum to sc^litSe orLinal 
reference spectra prior to subtracting the scaled spectL from tSe dJta. ^ 

IL ""f ^^^^"^ including the step of further scaling the spectra of 

co^eTESmedToThe'da^a.^^"^" ^ '^^'^ °^ ^ .odif ilatio^n ^d^^^^ 

?;ci:sLg-trst:p^of^^^^^^^ 

||f||g|C|. P.sub.kj spectra based on results of later sta'ges of'tt*?' "processing. 



22. The method as in claim 21 wherein said later stages of data Dro--ss^na and 
^"^^ys^s include principal components analysis (PCA) . pio..^ss^nq and 

23 The method as in claim 21 wherein said later stages of data processinq and 
^^■^^'/B^B include a technique of partial least squares (PLS) . ^ ^ """^ 

IL "method as in claim 1 or claim 2 or claim 3 or claim 4 or claim 17 includinq 
ZLI,? interactively displaying graphical output concerning which anaJytes aJd 
.hrf ?^ spectra, if any, are causing difficulties with respect to estimtion of 

Ire Ttlllllii'TZ interactive control over which reference spectra 

oraSL^ii^e^oJrStSf emJJ^yer'"'^' ^^^^'"^^"^ coefficients and wLghing 

^r-^/r^ ^" ^^^^"^ 24 including the step of further scaling the spectra of 

27. The method as in claim 1 or claim 2 or claim 3 or claim 4 or claim 17 ;,l«n 

^t^^x^^i^^ — siS^k.. 
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28. Apparatus for correcting input spectral data (X.sub.ki) derived from a sensor 
during a measurement, at least as to multiplicative errors, said apparatus 
comprising: 

jora/i^^^/^T^"^?''''!?^ ^ ^'3"^' representing spectral data subject to correction 
tor at least multiplicative errors; 

means for supplying a signal representing a first and primary reference spectrum 
(P.sub.ko) as a predetermined standard for such data; eierence spectrum 

Tp^sub^kl o^^P^suS\j)f^"^' representing at least one second reference spectrum 

means for estimating coefficients for a selected model for application to said incut 
spectral data, said input spectral data signal and fi^^ and sfJond reference 
spectrum signals being supplied to said estimating means; and 

e^MmJ^H^noffF^ ^° T^'^ estimating means for correcting said spectral data based on 
estimated coefficients at least as to multiplicative errors for producing a signal 
representing a linear additive structure for use in calibration. validaJf on anl 
determination by linear multivariate analysis . ■ ^^"-^^ 

30. Apparatus as in claim 28 wherein said means for supplying said signal 
representing at least one second reference spectrum supplies a signal which 
represents at least one spectrum representing expected influence of various 
undesired interferences (P.sub. jk) . 

31. Apparatus as in claim 28 wherein said means for supplying said signal 
representing at least one second reference spectrum supplies a signal which 
represents at least one spectrum representing expected influence of analytes of 
'"r^hf^/ V^""^ ^^"-^ ^"'^ including means for supplying a signal of at least 
one third reference spectrum representing expected influence of various undesired 
Jlggrgcji^^^said third reference spectrum signal also being supplied to said^' 

33^ Apparatus as in claim 28 or claim 29 or claim 30 or claim 31 wherein said 
estimating means estimates coefficients for a linear model . 

34 Apparatus as in claim 28 or claim 29 or claim 30 or claim 31 wherein said 
estimating means estimates coefficients for a non-linear model . 

35. In a system for analyzing a medium, said system having a spectrophotometric 
i^proJem^n^'co^^rj^ representing input spectral data (X.sub'ki,, the 

Spares: inclu'd^S-''"' '^'^ '^^^'^ multiplicative errors, said 

7^=nh^^''^^"^^^^^"^.^ ^^^"^^ representing a first and primary reference spectrum 
(P.sub.ko) as a predetermined standard for said spectral data; P «- um 

Tp^subS o^^P^sub^k?)f ^"^^ representing at least one second reference spectrum 

soJcJrlrdS' '^^?H^ r%"''''^f ^ selected model for application to said input 
spectral data, said input spectral signals and first and second reference spectral 
signals being supplied to said estimating means; and eterence spectral 

means responsive to said estimating means for correcting said spectral data based on 
estimated coefficients at least as to multiplicative errors for producing a siqnal 
representing a linear additive structure for use in calibration, ^al^datf on aS 
determination based on linear multivariate ana.l;jfsi3 . "«»i-ion ana 

36. The method for correcting input spectral data (X.sub.ki) derived from a 
measurement, by applying an additive model, and multiplicative model for the data in 
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a sequential fashion, cottprising the steps of: 

1) obtaining a set of spectral data Z from original input spectral data corrected 
for multiplicative effects by using a standard multiplicative scatter correction 
technique; 

2) estimating an additive model which takes into account spectral variations of 
analytes and interferences ; 

3) reconstructing spectral data Y without the estimated additive effects; 

4) estimating the multiplicative effects on Y; 

5) constructing a new matrix of corrected spectra Z from X; and 
repeating steps 2, 3 and 4 until convergence occurs. 

37^ The method of claim 36 wherein the additive effects are obtained from a 
different model than the multiplicative effects, 

38. The method of claim 37 wherein one model is a physical model and t he other model 
is a chemical model . = 

39. The method of claim 36 wherein step 4 includes the steps of: 
providing a first and primary reference spectrum; 

providing at least one second reference spectrum; 

estimating coefficients for a selected appropriate model to be applied to the input 
data based on said first and second reference spectra; and 

correcting said spectral data Y based on said estimated coefficients. 

40. The method for correcting input spectral data (X.sub.ki) derived from a 
measurement or for correcting the scores obtained by bilinear modeling of such data 
said method comprising the steps of: ' 

examining the measured spectral data as a combined additive and multiplicative 
structure such that 

X=T*P'+E 

where X=(x.sub.ik) is the matrix of spectral ordinates for i=l,2 . . , N objects, 
k=l,2 . . . , K wavelengths, T=(t. sub.il) is the matrix of scores for objects i,' 
1-1,2 . . . L representing bilinear factors Obtained from a bilinear mc^pl 
P=(P.sub.ki) are the loadings for objects i on bilinear factors 1, and E=(e sub ik) 
are the residuals between data X and model T*P'; 

decomposing the loadings into a function of a reference spectral and, optionally, a 
matrix of spectral components for analyte and interference phenomena; 

obtaining the additive and multiplicative effects for the input spectra from the 
coefficients from the loading decomposition and scores; and 

correcting said input spectra based on said obtained additive and multiplicative 

41. The method of claim 40 including the step of correcting said input spectral data 
simultaneously for inter£srent contributions and additive and multiplicative 
effects. 

42. The method for correcting input spectral data (x.sub.ki) derived from a 
measurement or for correcting the scores obtained by bilinear modelinq of such data, 
said method comprising the steps of: • ucua. 
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:??:ctuS sScVtSar^^ ^^^^ ^ ^""^'"^'^ ^-^^^^^^ multiplicative 

X=T*P'+E 

where X=(x.8ub^ik) is the matrix of spectral ordinates for i=l,2 . . . n objects 
k=l,2 . . K wavelengths, T=(t. sub.il) is the matrix of sco;es for objSs i 

I W-^ '^^P^^^^^'^i^g bilinear factors obtained from a bilinear model 
Iri^Z ^ ^^■/''^ loadings for objects i on bilinear factors 1, and E^T^b ik) 
are the residuals between data X and model T*P' ; le.suo.iKj 

decomposing the loadings into a function of a reference spectrum and, optionally a 
matrix of spectral components for analyte and interference phenomena?^' 

anSISrico^Llf'air '"''^'^ '^"^ ^ — ^-^-9- 

applying the obtained additive effects to the scores; 

e^i^i"^"^ multiplicative effect from the coefficient d from the reference 

obtaining scores corrected for additive and multiplicative effects using the 
multiplicative effect, or scores T and, optionally, the additive effect 

45. The method of claim 44 wherein the prediction procedure, using corrected 
spectral data, includes the steps of: ^ correctea 

using an independent set of data and applying factor loadings P to find a new set of 
spectral scores T; and 

determining the additive and multiplicative effects from the new spectral scores and 
o^tL^n^?^/ quantities of additive and multiplicative effects aSd those portions 
of the mterfei-ents and analytes which had been used in a prior calibration to 
obtain a corrected spectra; and ai-j-u" to 

using the corrected data and the calibration model coefficients from the linear 
model coefficients from the linear model to p^idlFt properties of interest 

IL of claim 44 wherein the prediction procedure, using corrected 

spectral scores, includes the steps of: 

Zeclrll icor^r^ranf' °' "^"'^'"^ '"^'^^'^^^ ' ^ of 

determining an S.sup.-l factor related to a matrix containing T*d for multiplicative 

jJlJtina S'aS'."'" ""T""^"^ "^""^"^ ' "^^"9 ^ " ^^^^ calibratioS 

relating to additive and mterferent effects to find 



W=S.sup.-l *T*U; where W is the matrix of the additive, multiplicative and 
S;isf^"' corrected scores which can be used as regressors ^in additive mixture 

Sctlrope^^Ss^'of J^tJJest^''''^''"" ^ coefficients from the linear model to 

47. Apparatus for correcting an input spectral data (X.sub.ki) signal derived from a 

::S":"'!nd"'' 1?"?'""'"'^ '° multiplicative errors, by the fit?Jng ofanlddl^rv: 
^!?.L and multiplicative model in sequential fashion, comprising: aaaitive 

orff^'.r?^"^ ^^'^ obtaining a signal representing a set of spectral data Z from the 
standard ™T fP^^tral data signal corrected for multiplicative effect by Jsi^g a 
standard multiplicative scatter correction technique; 
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second means for providing a signal representing estimate of an additive model which 
takes into account spectral variations of analytes and interferences ; 

third means for providing a signal representing the reconstructing of spectral data 
Y without the estimated additive effects; 

fourth means for providing a signal representing an estimate of the multiplicative 
effects on Y; 

fifth means for providing a signal representing the construction of a new matrix of 
corrected spectra Z from X; 

means for respectively providing said signal from said fifth means to said second, 
third and fourth means; and 

means responsive to the output signals of said second, third and fourth means for 
determining when convergence occurs. 

48. Apparatus for correcting an input spectral data (X.sub.ki) signal derived from a 
measurement or for correcting a signal representing the scores obtained by bilinear 
modeling of such data, comprising: 

means for providing a signal responsive to said input signal representing the 
measured spectral data as a combined additive and multiplicative structure such that 



X=T*P'+E 

where X=(x.sub.ik) is the matrix of spectral ordinates for i=l,z . . . N objects, 
k-1,2 . . . , k wavelengths, T=(t. sub.il) is the matrix of scores for objects i, 
1=1,2 . . . L representing bilinear factors obtained from a bilinear model , 
P=(P.sub.ki) are the loadings for objects i on bilinear factors 1 and E=(e.sub.ik) 
are the residuals between data X and model T*P» ; 

means for providing a signal representing decomposed loadings which have been 
decomposed into a function of a reference spectra and, optionally, a matrix of 
spectral components for analyte and interference phenomena; 

means for providing a signal representing the additive and multiplicative effects 
for the input spectra obtained from the coefficients from the loading decomposition 
and scores; and 

means for providing a signal representing the correction of said input spectra from 
said signal representing the additive and multiplicative effects. 

49. Apparatus for correcting an input spectral data (X.sub.ki) signal derived from a 
measurement or for correcting a signal representing the scores obtained by linear 
modeling of such data, comprising: 

means for providing a signal responsive to said input signal representing the 
measured spectral data as a combined additive and multiplicative structure such that 



X=T*P'+E 

where X=(x.sub.ik) is the matrix of spectral ordinates for i=l,2 . . . N objects, 
k=l,2 . . . , k wavelengths, T=(t. sub.il) is the matrix of scores for objects i, 
1=1,2 . . . L representing bilinear factors obtained from a bilinear "^oael, 
P=(P.sub.ki) are the loadings for objects i on bilinear factors 1 and E=Tersub. ik) 
are the residuals between data X and model T*P' ; 

means for providing a signal representing decomposed loadings which have been 
decomposed into a function of a reference spectra and, optionally, a matrix of 
spectral conponents for analyte and i?}t.?.3^^-erence phenomena; 
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means for providing a signal representing the additive effects for the scores T 
obtained from the offset corrected loadings and the scores T; 



means responsive to said signal representing the additive effects for providinq a 
llZt ^^P'^^^^"'^^"^ '^h^ application of the additive effects to the SJuiSa^iate 

means for providing a signal representing the multiplicative effect from the 
coefficient d from the reference spectrum; and ecrect trom the 

means for providing a signal representing scores corrected for additive and 

tSi addJtiJnflfct"'' "'"^ ^""^ nnaltiplicative effect, or scores T and, optionally 
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PUBLICATION-DATE: September 18, 2003 



File: PGPB 



INVENTOR- INFORMATION : 

NAME 

Ni, Jian 

Ruben, Steven M. 



CITY 

Germantown 
Brookeville 



STATE 

MD 

MD 



COUNTRY 

US 

US 



Sep 18, 2003 



RULE-47 



US-CL-CURRENT: 435/6; 435/ 320.1 , 435/325, 435/69.1, 530/350, 536/23. 2 

mm 




lof6 



10/23/2003 1:36 PM 



Record List Display 



httpy/westbrs:8()02/biD/gate.exe?^TOC&s..j^ 



O 3. Document ID: US 200301 14371 Al 

L32: Entry 3 of 11 File: PGPB 



Jun 19, 2003 



PG PUB -DOCUMENT -NUMBER; 20030114371 
PGPUB-FILING-TYPE: new 

DOCUMENT- IDENTIFIER: US 20030114371 Al 

JJbetaMr^^''''^^^^'''^'^^ encoding a novel human potassium channel beta-subunit , 

PUBLICATION-DATE: June 19, 2003 

INVENTOR- INFORMATION : 
NAME 

Feder, John N. 
Lee, Liana 
Chen, Jian 
Jackson, Donald 
Ramanathan, Chandra S. 
Siemers, Nathan O. 
Chang, Han 
Ryseck, Rolf -Peter 
Watson, Andrew J. 
Carroll, Pamela 



CITY 


STATE 


COUNTRY 


Belle Mead 


NJ 


US 


North Brunswick 


NJ 


US 


Princeton 


NJ 


US 


Lawrenceville 


NJ 


us 


Wallingford 


CT 


us 


Pennington 


NJ 


us 


Princeton Junction 


NJ 


us 


Ewing 


NJ 


us 


West Windsor 


NJ 


us 


Princeton 


NJ 


us 



RULE-47 



US -CL- CURRENT: 



: 514/12; 435/320.1, 435/325, 435/69.1, 530/350, 536/23,5 
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□ 4. Document ID: US 20030104487 Al 

L32: Entry 4 of 11 pil^. p^p^ 

PGPUB-DOCUMENT-NUMBER: 20030104487 
PGPUB-FILING-TYPE: new 

DOCUMENT- IDENTIFIER: US 20030104487 Al 

TITLE: Neuropeptide receptor and uses thereof 

PUBLICATION-DATE: June 5, 2003 

I NVENTOR - I NFORMAT I ON : 
NAME 

Fidock, Mark D. 
Robas, Nicola M. 



CITY 

Sandwich 
Sandwich 



STATE 



COUNTRY 
GB 
GB 



Jun 5, 2003 



RULE-47 



CL-CURRENT: 43S/7.2; 424/94 


.64 
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O 5. Document ID: US 20030036505 A 1 

L32: Entry 5 of 11 File: PGPB 

PGPUB-DOCUMENT-NUMBER: 20030036505 
PGPUB-FILING-TYPE: new 

DOCUMENT- IDENTIFIER: US 20030036505 Al 

TITLE: Signal transduction pathway component polynucleotides, polypeptides 
antibodies and methods based thereon P^xypepciaes, 



Feb 20, 2003 



PUBLICATION-DATE: February 20, 


2003 






INVENTOR- INFORMATION : 








NAME 


CITY 


STATE 


COUNTRY 


Barash, Steven C. 


Rockville 


MD 


US 


Ni, Jian 


Germantown 


MD 


US 


Ruben, Steven M. 


Olney 


MD 


US 


Rosen, Craig A. 


Laytonsville 


MD 


us 


Young, Paul E. 


Berkeley 


CA 


us 


Rohrschneider, Larry R. 


Seattle 


WA 


us 
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US-CL-CURRENT: 514/12; 435/320^, 435/325, 435/6, 435/69.1, 530/350, 536/23.5 
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Feb 13, 2003 



O 6. Document ID: US 20030032786 Al 

L32: Entry 6 of 11 File: pgpb 

PGPUB-DOCUMENT-NUMBER: 2003 0032786 
PGPUB-FILING-TYPE: new 

DOCUMENT- IDENTIFIER: US 20030032786 Al 

l^l^f'-J^^^y^^^^^otide encoding a novel human potassium channel beta-subunit , 
K+DetaM2 ' 



PUBLICATION-DATE: February 13, 


2003 






I NVENTOR - I NFORMAT I ON : 








NAME 


CITY 


STATE 


COUNTRY 


Chang, Han 


Princeton Junction 


NY 


US 


Chen, Jian 


Princeton 


NJ 


US 


Feder, John 


Belle Mead 


NJ 


US 


Jackson, Donald 


Lawrenceville 


NJ 


US 


Lee, Liana 


North Brunswick 


NJ 


US 


Ramanathan, Chandra S. 


Wallingford 


CT 


US 


Siemers, Nathan 0. 


Pennington 


NJ 


US 


Carroll, Pamela 


Princeton 


NJ 


us 
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US -CL- CURRENT: 536/23.1 
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□ 7. Document ID: US 20020187950 Al 



L32: Entry 7 of 11 

PGPUB-DOCUMENT-NUMBER: 20020187950 
PGPUB-FILING-TYPE: new 

DOCUMENT-IDENTIFIER: US 20020187950 Al 
TITLE: Keratinocyte derived interfe ron 
PUBLICATION-DATE: December 12, 2002 



INVENTOR- INFORMATION : 
NAME 

LaFleur, David W. 
Moore, Paul A. 
Ruben, Steven M. 



CITY 

Washington 
Germantown 
Olney 



File: PGPB 



STATE 
DC 
MD 
MD 



COUNTRY 

US 

US 

us 



US -CL- CURRENT 



Dec 12, 2002 



RULE- 4 7 



: 514/44; 424/85^, 435/320J,. 435/325, 435/69.51. 530/351, _536/23.5 



i Title i cfatio» 1 F,ont j R€^;;;;i"cSs;fi;;3^o»T Ti.f...... 



equenees Attachment 



Nov 14, 2002 



O 8. Document ID: US 20020169292 A 1 

L32: Entry 8 of 11 Pil^^ p^.p3 

PGPUB-DOCUMENT-NUMBER: 20020169292 
PGPUB-FILING-TYPE: new 

DOCUMENT- IDENTIFIER: US 20020169292 Al 

TITLE: Cystine knot growth factor mutants 

PUBLICATION-DATE: November 14, 2002 

INVENTOR- INFORMATION : 
NAME 

„ . ^ ^ ^^'^^ STATE COUNTRY RULE-47 

Weintraub, Bruce D. Rockville md US 

Szkudlinski, Mariusz W. Potomac md us 

US-CL-CURRENT: 530/397; .435/ 320.1 , 435/325, .135/69,4, 536/23.5 



O 9. Document ID: US 6472512 Bl 

L32: Entry 9 of 11 pil ^SPT 

tixe. UbPT Oct 29, 2002 

US-PAT-NO: 6472512 

DOCUMENT-IDENTIFIER: US 6472512 Bl 
TITLE: Keratinocyte derived interferon 
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.PCPBJPAB^PABJDWPI,TDBD&ESNAME^ 



DATE- ISSUED: October 29, 2002 

I NVENTOR - I NFORMAT I ON : 
NAME 

LaFleur; David W. 
Moore; Paul A. 
Ruben; Steven M. 



CITY 

Washington 
Germantovm 
Olney 



STATE 
DC 
MD 
MD 



ZIP CODE 



COUNTRY 



s%?39fT^'' 435/335, 43S/7.92. 530/ 388.15 , 530/3S9.2, 





^::|:C/tatlor 


> j Ffont ] 










Sequences | Attachments i 


pDraw. Deso j 

















O 10. Document ID: US 6433145 Bl 

L32: Entry 10 of 11 
US-PAT-NO: 6433145 

DOCUMENT- IDENTIFIER: US 6433145 Bl 

** See image for Certificate of Correction ** 

TITLE: Keratinocyte derived interferon 
DATE- ISSUED: August 13, 2002 



File: USPT 



Aug 13, 2002 



INVENTOR- INFORMATION : 
NAME 

LaFleur; David W. 
Moore; Paul A. 
Ruben; Steven M. 



CITY 

Washington 
German town 
Olney 



STATE 
DC 
MD 
MD 



US-CL-CURRENT: 530/351; 424 / 85.4 , 435/7.1, 530/350 



ZIP CODE 



COUNTRY 



1 "fu8 [ fitie 1 Citattoi 


vr.FfCint j Rsvt6iAi [ Claj^'if 


io3tion { Date ] R 


efefsnoe ]■ Secjuences 1 
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O 11. Document ID: US 6375634 Bl 

L32: Entry 11 of 11 File: USPT 

US-PAT-NO: 6375634 

DOCUMENT- IDENTIFIER: US 6375634 Bl 



Apr 23, 2002 



TITLE: APPARATUS AND METHOD TO ENCAPSULATE, KILL AND REMOVE MALIGNANCIES INCLUDING 
SELECTIVELY INCREASING ABSORPTION OF X-RAYS AND INCREASING FrS-rSiS^ ^^^^ 
RESIDUAL TUMORS TARGETED BY IONIZING AND NON-IONIZING RADIATION THE^ 

DATE-ISSUED: April 23, 2002 



INVENTOR- INFORMATION : 
NAME 

Carroll; Robert G. 



CITY 
Largo 



STATE 
FL 



ZIP CODE 



COUNTRY 
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jCYROSTAT-CUT 
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■1 oil 
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II "I 




Search: l fj | Refine Search 



.El 



Search History 
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Set Name Query Hit Count Set Name 

side by side result set 

DB=USPT.PGPB.JPAB.EPAB.DWPI.TDBD; PLUR=YES; OP=ADJ 

L4 L3 and (cyrostat$6) 0 L4 
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Term 
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Oil 
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1 2345j 


NODA 


j 3025111 


NODAB 




NODABC 




NODABCD 








NODABCDIJ 






i 2 


ilNODABCDIJL 






\ 1 


JNODABCIJ 






i^ODABILITY 
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r ::zzi: 


i|NODABLE 
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Set Name Query 

side by side 

DB=USPT,PGPB,JPAB.EPAB.DWPI.TDBD: PLUR=YES; OP=ADJ 
LIO L9 and (nod$7) 

L8 and (instab$8 or unstabl$4 or inhomogene$8 or non-uniform$6 or 
L9 nonunifonn$6 or distrub$6 or disrupt$6 or flux$2 or fluct$9 or 
vari$9) 

J g L7 and (corq)ut$9 or proces$8 or analyz$6 or analys$4 or calculat$8 
— or model$6 or simulat$9) 

L7 L6 and (eddy or vfl3rat$9 or ghost$4 or blur$6 or alias$6 or stiff$6) 
L6 L5 and (cryo$9) 
L5 L4 and (damp$6) 

L4 L3 and (motion$4 or movement or mov$4) 
L3 L2 and (gradient) 

L2 LI and (temperature or thennal$4 or heat$4) 
LI ((magnetic adj resonance) or MRI or NMR) 



Hit Count Set Name 
result set 



13 LIO 



45 L9 



45 

45 
58 
348 
7221 
24182 
126088 
158767 



L8 

LZ 
L6 
L5 
L4 
L3 
L2 
LI 
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Term 

||(30i^^ 

|(L30NOfL31).USP^^^^^ 



Documents i| 

~ Toll 



Database: 




Search: 



I L34 



Refine Search 



Search History 



DATE: Thursday, October 23, 2003 Printable Copy Create Case 



Set Name Query 
side by side 

DB=USPT,PGPB.JPAB.EPAB,DWPI,TDBD: PLUR=YES; OP=ADJ 
L34 LSOnotBl 
L31 not 132 
L31 and(nod$7) 

L30 and (motion$4 or movement or mov$4) 

L29 and ((analysis or analyz$6 or process$4 or comput$9) with 
(vftratS? or stress$4 or strain$5 or danp$7 or stiff$5) ) 

L28 and ((model$4 or simulat$5) with (instability or unstable or 
inhomogeneit$4 or eddy or disturbs? or interfer$6 or non-uniform$5 
or non-uniform$5)) 

L28 L27 and (cool$4 or cyro$9) 



Hit Count Set Name 



133 
L32 
L31 

L30 
L29 





result set 


10 


[.34 


21 


1.33 


11 


L32 


32 


L31 


42 


1.30 



161 L29 



6066 L28 
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L27 


L26 and (instability or unstable or inhomogeneit$4 or eddy or 
disturb$7 or interfer$6 or non-uniform$5 or non-unifonn$5) 




i 0 7 


126 


L25 and (analysis or analyz$6 or process$4) 


16873 


L26 


L25 


LI 2 and (vibrat$7 or stress or strain or danp$7 or stiff$5) 


17216 


L25 


L24 


L23 and L20 


9 


L24 




{((^Odf^OO l^94Ani I'lOA/'Xno llOA/^n^ \inA/inA iio/i/iAc \^^A/'-if\^ 
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1324/321 1324/322 )!.CCLS.)) 


6736 


L23 


L22 


L21 and (motion or movement or mov$3) 


19 


L22 


in 


L20 and ((model$4 or simulat$5) with (instability or unstable or 
inhomogeneit$4 or eddy or disturb$7 or interfer$6)) 


22 


L21 


T OA 

LZV 


LI 9 and ((analysis or analyz$6 or process$4) with (vibrat$7 or stress 
or strain or damp$7 or stif&ess) ) 


402 


L20 


L19 


LIS and ((analysis or analyz$6 or process$4) with (instability or 
unstable or inhomogeneit$4 or eddy or disturb$7 or interfer$6)) 




T 1 O 

L.iy 


LI 8 


LI 6 and (cool$4 or cyro$9) 


3521 


L18 


L17 


LI 6 and (cyrostat) 


1 


LI 7 


L16 


LI 5 and (transfer$4) 


6292 


L16 


L15 


LI 4 and (instability or unstable or inhomogeneit$4 or eddy or 
disturb$7 or interfer$6) 






Li4 


LI 3 and (analysis or analyz$6 or process$4) 


13050 


L14 


LI 3 


LI 2 and (vibrat$7 or strss or strain or damp$7 or stiffiiess) 


13293 


LI 3 


L12 


LI 1 and (model$4 or simulat$5) 


42595 


Ul 


Lll 


((magnetic adj resonance) or MRI or NMR) 


158767 


Lll 


LIO 


L9 and (nod$7) 


13 


LIO 


L9 


L8 and (instab$8 or unstabl$4 or inhomogene$8 or non-uniform$6 or 
nonuniform$6 or distrub$6 or disrupt$6 or flux$2 or fluct$9 or 
vari$9) 


45 


13 


L8 


L7 and (comput$9 or proces$8 or analyz$6 or analys$4 or calculat$8 
or model$6 or simulat$9) 




■Lo 


U 


L6 and (eddy or vibrat$9 or ghost$4 or blur$6 or alias$6 or stiff$6) 


45 


L7 


L6 


L5 and (cryo$9) 


58 


L6 


L5 


L4 and (damp$6) 


348 


L5 


L4 


L3 and (motion$4 or movement or mov$4) 


7221 


L4 


L3 


L2 and (gradient) 


24182 


L3 


L2 


LI and (temperature or thermal$4 or heat$4) 


126088 


Ul 


LI 


((magnetic adj resonance) or MRI or NMR) 


158767 


Li 
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Search Results - Record(s) 1 through 13 of 13 returned. 



□ 1. Document ID: US 20030166277 Al 
LIO: Entry 1 of 13 

PGPUB-DOCUMENT-NUMBER: 20030166277 
PGPUB-FILING-TYPE: new 

DOCUMENT -IDENTIFIER: US 20030166277 Al 
TITLE: Targeted vaccine delivery systems 
PUBLICATION-DATE: September 4, 2003 



INVENTOR- INFORMATION : 
NAME 

Zauderer, Maurice 
Smith, Ernest S. 



CITY 

Pittsford 
Rochester 



File: PGPB 



STATE 

NY 

NY 



US -CL- CURRENT: 435/372; 424 / 178.1 , 530 / 391.: 



COUNTRY 

US 

US 



Sep 4, 2003 



RULE-47 



FuB j Title I Citation i 



O 2. Document ID: US 20030157113 Al 

LIO: Entry 2 of 13 File: PGPB 

PGPUB-DOCUMENT-NUMBER: 20030157113 
PGPUB-FILING-TYPE: new 

DOCUMENT- IDENTIFIER: US 20030157113 Al 

TITLE: Compositions and methods for treatment of neoplastic disease 
PUBLICATION-DATE: August 21, 2003 
INVENTOR- INFORMATION : 



Aug 21, 2003 



NAME 


CITY 


STATE 


COUNTRY RULE-47 


Terman, David S. 


Pebble Beach 


CA 


US 


US -CL- CURRENT: 424/184 . 1 


; 435/346 






FuB j Title i Citation j Front :J 











□ 3. Document ID: US 20030153063 A 1 

LIO: Entry 3 of 13 File: PGPB 



Aug 14, 2003 
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PGPUB-DOCUMENT-NUMBER: 20030153063 
PGPUB-FILING-TYPE: new 

DOCUMENT -IDENTIFIER: US 20030153063 Al 

TITLE: Novel human G-protein coupled receptor, HGPRBMYll, expressed highly in heart 
and variaat:s thereof 



PUBLICATION-DATE: August 14, 2003 
INVENTOR- INFORMATION : 



NAME 


CITY 


STATE 


COUNTRY 


Feder, John N. 


Belle Mead 


NJ 


US 


Nelson, Thomas C. 


Lawrenceville 


NJ 


US 


Ramanathan, Chandra S. 


Wallingford 


CT 


us 


Cacace, Angela M. 


Clinton 


CT 


us 


Barber, Lauren E. 


Griswood 


CT 


us 



RULE -47 



US-CL-CURRENT: 435/226; 435 / 320.1 , 435/325, 435/6, 435 / 69.1 , 536 / 23.2 



n 4. Document ID: US 20030107376 A 1 

LIO: Entry 4 of 13 File: PGPB Jun 12, 2003 

PGPUB-DOCUMENT-NUMBER: 20030107376 
PGPUB-FILING-TYPE: new 

DOCUMENT- IDENTIFIER: US 20030107376 Al 
TITLE: Magnetic resonance imaging apparatus 
PUBLICATION-DATE: June 12, 2003 
INVENTOR- INFORMATION : 

CITY STATE COUNTRY RULE-47 

Yasuhara, Yasutake Nasu-gun jp 

US-CL-CURRENT: 324/318; 324/320, 324/321 




D 5. Document ID: US 20030099608 A 1 

LIO: Entry 5 of 13 File: PGPB May 29, 2003 

PGPUB-DOCUMENT-NUMBER: 20030099608 
PGPUB-FILING-TYPE: new 

DOCUMENT- IDENTIFIER: US 20030099608 Al 
TITLE: Human cytokine receptor 
PUBLICATION-DATE: May 29, 2003 
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s..^|^PT, 



INVENTOR- INFORMATION : 

NAME 

Presnell, Scott R. 
Xu, Wenfeng 
Kindsvogel, Wayne 
Chen, Zhi 
Hughes, Steven D. 



LI 1 I 


STATE 


COUNTRY 


Tacoma 


WA 


US 


Mukilteo 


WA 


US 


Seattle 


WA 


US 


Bellevue 


WA 


US 


Seattle 


WA 


us 



RULE-47 



US-CL-CURRENT: 424 / 85.1 ; 43S / 320.1 , 435/325, 435/69.S, 530/351, 536/23.5 




O 6. Document ID: US 20030077706 Al 

LIO: Entry 6 of 13 File: PGPB 

PG PUB -DOCUMENT -NUMBER: 2003 0077706 
PGPUB-FILING-TYPE: new 

DOCUMENT- IDENTIFIER: US 20030077706 Al 
TITLE: Mouse cytokine receptor 
PUBLICATION-DATE: April 24, 2003 
I NVENTOR - I NFORMAT I ON : 

CITY STATE COUNTRY RULE-47 

Presnell, Scott R. 
Xu, Wenfeng 
Kindsvogel, Wayne 
Chen, Zhi 

US-CL-CURRENT: 115/69.1; 135/320.1 435/525, 135/6, 530/350, 536/23.5 



CITY 


STATE 


COUNTRY 


Tacoma 


WA 


US 


Mukilteo 


WA 


US 


Seattle 


WA 


US 


Bellevue 


WA 


US 



Apr 24, 2003 



O 7. Document ID: US 20030006773 Al 

LIO: Entry 7 of 13 File: PGPB Jan 9, 2003 

PGPUB-DOCUMENT-NUMBER: 2003 0006773 
PGPUB-FILING-TYPE: new 

DOCUMENT- IDENTIFIER: US 20030006773 Al 

TITLE: Magnetic resonance apparatus having a basic field magnet 
PUBLICATION-DATE: January 9, 2003 
INVENTOR- INFORMATION : 

CITY STATE COUNTRY RULE-47 

Ries, Guenter Erlangen DE 
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O 8. Document ID: US 20020177551 Al 

LIO: Entry 8 of 13 File: PGPB 

PGPUB-DOCUMENT-NUMBER: 20020177551 
PGPUB-FILING-TYPE: new 

DOCUMENT -IDENTIFIER: US 20020177551 Al 

TITLE: Compositions and methods for treatment of neoplastic disease 
PUBLICATION-DATE: November 28, 2002 



I NVENTOR - I NFORMAT I ON : 
NAME 

Terman, David S. 



CITY 

Pebble Beach 



STATE COUNTRY 
CA US 



Nov 28, 2002 



RULE-47 



US -CL- CURRENT: 514/12; 435/325, 53C7350 
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□ 9. Document ID: US 20010022515 Al 

LIO: Entry 9 of 13 File: PGPB 

PGPUB-DOCUMENT-NUMBER: 20010022515 
PGPUB-FILING-TYPE: new 

DOCUMENT- IDENTIFIER: US 20010022515 Al 

TITLE: Magnetic resonance imaging apparatus 

PUBLICATION-DATE: September 20, 2001 

INVENTOR- INFORMATION : 
NAME 

Yamashita, Masatoshi 
Ishii, Manabu 
Sakakura , Yoshi tomo 
Takamori , Hiromitsu 
Nakabayashi , Kazuto 
Hamamura, Yoshinori 
Mitsui, Shinji 
Yasuhara , Yasu take 



CITY 


STATE 


COUNTRY 


Otawara-shi 




JP 


Otawara-shi 




JP 


Nasu-gun 




JP 


Otawara-shi 




JP 


Nasu-gun 




JP 


Otawara-shi 
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